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' This paper investigates joint wireless information and energy transfer in a two-user MIMO interfer- 

ence channel, in which each receiver either decodes the incoming information data (information decoding, 
ID) or harvests the RF energy (energy harvesting, EH) to operate with a potentially perpetual energy 
supply. In the two-user interference channel, we have four different scenarios according to the receiver 
mode - (IDi, ID 2 ), {EH X , EH 2 ), (EHi, ID 2 ), and {ID X , EH 2 ). While the maximum information 
O ' bit rate is unknown and finding the optimal transmission strategy is still open for (ID X , ID 2 ), we have 

derived the optimal transmission strategy achieving the maximum harvested energy for {EH\, EH 2 ). For 

CN . 

(EHi, ID 2 ), and (IDi, EH 2 ), we find a necessary condition of the optimal transmission strategy and, 

CO ■ 

Q\ • accordingly, identify the achievable rate-energy (R-E) tradeoff region for two transmission strategies 



X 



that satisfy the necessary condition - maximum energy beamforming (MEB) and minimum leakage 



CO . beamforming (MLB). Furthermore, a new transmission strategy satisfying the necessary condition - 

O ■ 

^f*) . signal-to-leakage-and-energy ratio (SLER) maximization beamforming - is proposed and shown to exhibit 

a better R-E region than the MEB and the MLB strategies. Finally, we propose a mode scheduling method 
to switch between (EHi, ID 2 ) and (IDi, EH 2 ) based on the SLER. 
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I. Introduction 

Over the last decade, there has been a lot of interest to transfer energy wirelessly and recently, radio- 
frequency (RF) radiation has become a viable source for energy harvesting. It is nowadays possible to 
transfer the energy wirelessly with a reasonable efficiency over small distances and, furthermore, the 
wireless sensor network (WSN) in which the sensors are capable of harvesting RF energy to power their 
own transmissions has been introduced in industry ( [l]-[4] and references therein). 

The energy harvesting function can be exploited in either transmit side [5]-[9] or receive side [10]- 
[13]. For the energy harvesting transmitter, energy harvesting scheduling and transmit power allocation 
have been considered and, for the energy harvesting receiver, the management of information decoding 
and energy harvesting has been developed. Furthermore, because RF signals carry information as well as 
energy, "joint wireless information and energy transfer" in conjunction with the energy harvesting receiver 
has been investigated [10]— [13]. That is, previous works have studied the fundamental performance limits 
and the optimal transmission strategies of the joint wireless information and energy transfer in the cellular 
downlink system with a single base station (BS) and multiple mobile stations (MSs) [12] and in the 
cooperative relay system [13] and in the broadcasting system [10], [11] with a single energy receiver and 
a single information receiver when they are separately located or co-located. 

There have been very few studies of joint wireless information and energy transfer on the interference 
channel (IFC) models [14]— [16]. In [14], [15], the authors have considered a two-user single-input single- 
output (SISO) IFC and derived the optimal power scheduling at the energy harvesting transmitters that 
maximizes the sum-rate given harvested energy constraints. In [16], the authors have investigated joint 
information and energy transfer in multi-cell cellular networks with single-antenna BSs and single-antenna 
MSs. To the best of the authors' knowledge, the general setup of multiple-input multiple-output (MIMO) 
IFC models accounting for joint wireless information and energy transfer has not been addressed so far. 

As an initial step, in this paper, we investigate a joint wireless information and energy transfer in a 
two-user MIMO IFC, where each receiver either decodes the incoming information data (information 
decoding, ID) or harvests the RF energy (energy harvesting, EH) to operate with a potentially perpetual 
energy supply. Because practical circuits and hardware that harvest energy from the received RF signal 
are not yet able to decode the information carried through the same RF signal [10], [11], [17], we assume 
that the receiver cannot decode the information and simultaneously harvest energy. It is also assumed 
that the two (Tx l,Tx 2) transmitters have knowledge of their local CSI only, i.e. the CSI corresponding 
to the links between a transmitter and all receivers (Rx 1, Rx 2). In addition, the transmitters do not 
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share the information data to be transmitted and their CSI and, furthermore, the interference is assumed 
not decodable at the receiver nodes as in [18]. That is, Tx 1 (Tx 2) cannot transfer the information 
to Rx 2 (Rx 1). In a two-user IFC, we then have four different scenarios according to the Rx mode 
- (IDx, ID 2 ), (EH U EH 2 ), (EH i, ID 2 ), and (IDx, EH 2 ). Because, for (ID U ID 2 ), the maximum 
information bit rate is unknown and finding the optimal transmission strategy is still an open problem 
in general, we investigate the achievable rate when a well-known iterative water-filling algorithm [19]- 
[21] is adopted for (IDx, ID 2 ) with no CSI sharing between two transmitters. For (EH\, EH 2 ), we 
derive the optimal transmission strategy achieving the maximum harvested energy. Because the receivers 
operate in a single mode such as (ID\, ID 2 ) and (EH\, EH 2 ), when the information is transferred, no 
energy is harvested from RF signals and vice versa. For (EH\, ID 2 ) and (ID\, EH 2 ), the achievable 
energy-rate (R-E) trade-off region is not easily identified and the optimal transmission strategy is still 
unknown. However, in this paper, we find a necessary condition of the optimal transmission strategy, 
in which one of the transmitters should take a rank-one energy beamforming strategy with a proper 
power control. Accordingly, the achievable R-E tradeoff region is identified for two different rank-one 
beamforming strategies - maximum energy beamforming (MEB) and minimum leakage beamforming 
(MLB). Furthermore, we also propose a new transmission strategy that satisfies the necessary condition - 
signal-to-leakage-and-energy ratio (SLER) maximization beamforming. Note that the SLER maximizing 
approach is comparable to the signal-to-leakage-and-noise ratio (SLNR) maximization beamforming [22], 
[23] which has been developed for the multi-user MIMO data transmission, not considering the energy 
transfer. The simulation results demonstrate that the proposed SLER maximization strategy exhibits wider 
R-E region than the conventional transmission methods such as MLB, MEB, and SLNR beamforming. 
Finally, we propose a mode scheduling method to switch between (EH\, ID 2 ) and (ID\, EH 2 ) based 
on the SLER that further extends R-E tradeoff region. 

The rest of this paper is organized as follows. In Section |Il] we introduce the system model for 
two-user MIMO IFC. In Section [TTlJ we discuss the transmission strategy for two receivers on a single 
mode, i.e. (ID\, ID 2 ) and (EH\, EH 2 ). In Section |IVJ we derive the necessary condition for the 
optimal transmission strategy and investigate the achievable rate-energy (R-E) region for (EH\, ID 2 ) 
and (ID±, EH 2 ) and, in Section [V] propose the SLER maximization strategy. In Section [VTI and Section 
IVII1 we provide several discussion and simulation results, respectively, and in Section IVIIII we give our 
conclusions. 

Throughout the paper, matrices and vectors are represented by bold capital letters and bold lower-case 
letters, respectively. The notations (A)^, (A)t, (A)j, [A],, tr(A), and det(A) denote the conjugate 
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Fig. 1. Two-user MIMO IFC in (EHi, ID 2 ) mode. 



transpose, pseudo-inverse, the ith row, the ith column, the trace, and the determinant of a matrix A, 
respectively. The matrix norm ||A|| and ||A||_p denote the 2-norm and Frobenius norm of a matrix A, 
respectively, and the vector norm ||a|| denotes the 2-norm of a vector a. In addition, (a) + = max(a, 0) 
and A y means that a matrix A is positive semi-definite. Finally, 1m denotes the M X M identity 
matrix. 

II. System model 

We consider a two-user MIMO IFC system where two transmitters, each with M t antennas, are 
simultaneously transmitting their signals to two receivers, each with M r antennas, as shown in Fig. 
[TJ Note that each receiver can either decode the information or harvest energy from the received signal, 
but it cannot execute the information decoding and energy harvesting at the same time due to the hardware 
limitations. That is, each receiver can switch between ID mode and EH mode at each frame or time slotj^ 
We assume that the transmitters have perfect knowledge of the CSI of their associated links (i.e. the 
links between a transmitter and all receivers) but do not share those CSI between them. In addition, 
M t = M r = M for simplicity, but it can be extended to general antenna configurations. Assuming a 
frequency flat fading channel, which is static over several frames, the received signal yj G C xl f or 

'Note that the switching criterion between ID mode and EH mode depends on the receiver's condition such as the available 
energy in the storage and the required processing or circuit power. In this paper, we focus on the achievable rate and harvested 
energy obtained by the transferred signals from both transmitters in the IFC according to the different receiver modes. The mode 
switching policy based on the receiver's condition is left as a future work. We assume that the mode decided by the receiver is 
sent to both transmitters through the zero-delay and error- free feedback link at the beginning of the frame. 



May 1, 2013 



DRAFT 



5 



i = 1, 2 can be written as 

yi = Hnxi + Hi 2 x 2 + ni, 

y 2 = H 2 ixi + H 22 x 2 + n 2 , (1) 

where rij G C xl is a complex white Gaussian noise vector with a covariance matrix Im and Hjj G 
C MxM is the normalized frequency-flat fading channel from the jth transmitter to the zth receiver such 
as Ejfc=i \ h ?j k) \ 2 = &ij M [ 24 ]- Here, hfj k) is the (l,k)Hx element of and o^- G [0, 1]. We assume 
that has a full rank. The vector Xj G C Mxl is the transmit signal, in which the independent messages 
can be conveyed, at the jth transmitter with a transmit power constraint for j = 1 and 2 as 

^[llxjH 2 ] < P for j = 1 and 2. (2) 

When the receiver operates in ID mode, the achievable rate at zth receiver, Ri, is given by [19] 

Ri = logdet(I M + HgRljHaQi), (3 ) 

where R_; indicates the covariance matrix of noise and interference at the ith receiver, i.e., 

1 = Ijw + Hi 2 Q 2 H^, 
R- 2 = lM + H 2 iQiHg. 

Here, Qj = i?[xjx|^] denotes the covariance matrix of the transmit signal at the jth transmitter and, 
from ©, tr(Qj) < P. 

For EH mode, it can be assumed that the total harvested power Ei at the ith receiver (more exactly, 
harvested energy normalized by the baseband symbol period) is given by 

Ei = tiE[\\ yi \\ 2 \ 

= Qtr ^H^Q.H^ + ImV (4) 

where Q denotes the efficiency constant for converting the harvested energy to electrical energy to be 
stored [3], [10]. For simplicity, it is assumed that d = 1 and the noise power is negligible compared to 
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the transferred energy from each transmitters j^| That is, 

= tr (HuQiHg) + tr (H i2 Q 2 Hg) 

= En + E^, (5) 

where Eij = tr ^HyQjH^ denoting the energy transferred from the jth transmitter to the ith receiver. 

Interestingly, when the receiver decodes the information data from the associated transmitter under the 
assumption that the signal from the other transmitter is not decodable [18], the signal from the other 
transmitter becomes an interference to be defeated. In contrast, when the receiver harvests the energy, it 
becomes a useful energy-transferring source. Fig. Q] illustrates an example of the receiving mode (EHi, 
ID2), where the interference 1 (dashed red line) should be reduced for ID, while the interference2 (dashed 
green line) be maximized for EH. In what follows, for four possible receiving modes, we investigate the 
achievable rate-harvested energy tradeoff. In addition, the corresponding transmission strategy (more 
specifically, transmit signal design) is presented. 

III. TWO RECEIVERS ON A SINGLE MODE 

A. Two IDs: maximum achievable sum rate 

For the scenario (IDi, ID 2 ), it is desirable to obtain the maximum achievable sum rate. That is, the 
problem can be formulated as follows: 



(PI) maximize Y2i=x &i (6) 

subject to tr(Qj) < P, Qj >z for j = 1, 2, (7) 

The solution of (PI) has been extensively considered in many previous communication researches 
[19]— [21], where the iterative water-filling algorithms have been developed to maximize the achievable 
rate in a distributed manner with no CSI sharing between the transmitters. This is briefly summarized in 
Algorithm 1: 

Algorithm 1. Iterative Water-filling: 

2 In this paper, we assume the system operates in the high signal-to-noise ratio (SNR) regime, which is also consistent with 
the practical wireless energy transfer requires a high-power transmission, but we also discuss the low SNR regime in Section 
IVilas well. 
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1) Initialize n = and G Qp for j = 1, 2, where 



Q P = {Q G C MxM : Q >r 0, ir(Q) = P}. (8) 

,1 



2) For n = : N rnax , where N max is the maximum number of iterations 



Update Q^ n+1 ^ for j = 1, 2 as follows: 



^(Hjj , R^j , P) , if R£J is updated, 
Q^™' ) , otherwise, 

where R_^. indicates the covariance matrix of noise and interference in the jth receiver at the nth 
iteration, i.e., 

RL1=lAf + H 12 Q^Hf 2 , 
R^?2 = Im + H2iQ^H^. 

(n) 

Note that Ri- is measured at each receiver similarly to the way it has been done in [19] and, 

(n) 

furthermore, is computed at the receiver and reported to the transmitter through the zero-delay 
and error-free feedback link. 
3) Finally, = Q^ m " +1 for j = 1, 2. 

Here, WF{) denotes the water-filling operator given as [19]: 

WF(H ih R, P) = Vi(fHlM ~ D^J+Uf , (10) 

where Uj and Dj are obtained from the eigenvalue decomposition of H^R _1 Hjj. That is, H^R _1 Hjj = 
UjDjU^, and //j denotes the water level that satisfies the transmit power constraint as tr{{^ii.M — 
D f ri)+} = P. 

In the scenario (ID\, ID2), because both receivers decode the information, the harvested energy 
becomes zero. 

B. Two EHs: maximum harvested energy 

For the scenario {EH\, EH2), both receivers want to achieve the maximum harvested energy. That is, 
the problem can be formulated as: 

(P2) maximize Ei (11) 

subject to tr(Qj) < P, Qj h for j = 1, 2, (12) 

3 Generally, Nmax = 20 is sufficient for the solutions to converge. 
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The following proposition gives the optimal solution for the problem (P2). 

Proposition 1: The optimal Qj for (P2) has a rank equal to one and is given as Qj = P[Vj]i[Vj]^, 

.That is, Hj = UjEjVf , 



Hi, 



H 



2j 



where Vj is a M x M unitary matrix obtained from the S VD of Hj = 

where £j = diag{&j t i, &j t M } with ct^i > ... > (7 j,Af ■ 
Proof: From d5), 

2 2/2 

i=l i=l \j=l 

j=i \i=i j 

2 

= E* r ( fi iQi S f) < 13 ) 

i=i 

Note that the covariance matrix Qj can be written as Qj = VjD 2 V^ where Vj is a Af x M unitary 
matrix and D| = diag{d 2 ^ ...,d 2 M } with d£ m < P. Because ir(AB) = tr(BA) for A G C mxn 

and B G C nxm , (|T3T > can be rewritten as 

2 2 2 Af 

£3» = X> (DfVf Hf H,Vj) = £ E 4J|Hi[Vj] m || 2 . (14) 
Because d 2 m < P, 

Af 

^4 m ||Hj[Vj] m || 2 <P max IIH^V.UH 2 . (15) 

* — ' m=l,...Af 
m=l 

Here, the equality holds when d? jm , = P for ml = arg max ||Hj[Vj] m || 2 and d 2 - = for m / m! , 

J ' m=l,...M J ' 

which implies that Qj has a rank equal to one and accordingly, it is given as Qj = P\V j\m> [Vj]^,. 
Note that 

||Hj[V,] m ,|| 2 <a] A , (16) 
where the equality holds when [Vj] m / = \Vj]i- Therefore, from ([T5l and (TTol , (PT4l is bounded as 

2 2 Af 

E^ = E E 4ml|H,[Vj] m || 2 < P(aj A +a 2 2A ), 

i=l j=l m=l 

and the equality holds when Qj = P[Vj]i[Vj]f . ■ 
Note that each transmitter can design the transmit covariance matrix Qj such that the transferred energy 
from each transmitter is maximized without considering other transmitter's channel information and 



May 1, 2013 



DRAFT 



9 



transmission strategy. That is, thanks to the energy conservation law, each transmitter transfers the energy 
through its links independently. 

From Proposition Q] the transmit signal on each transmitter can be designed as Xj = vPfVjjiSy, 
where Sj is any random signal with zero mean and unit variance. Because both receivers harvest the 
energy and are not able to decode the information, the achievable rate becomes zero. 

IV. One ID receiver and One EH receiver 

In this section, without loss of generality, we will consider (EH±, ID 2 ) - the first receiver harvests 
the energy and the second decodes information. The transmission strategy described below can also be 
applied to (ID\, EH 2 ) without difficulty. Note that energy harvesting and information transfer occur 
simultaneously in the IFC, and accordingly, the achievable rate-energy region is not trivial compared to 
the scenarios {EH\, EH 2 ) and (JDi, ID 2 ). 

A. A necessary condition for the optimal transmission strategy 

Because information decoding is done only at the second receiver, by letting R = R 2 and E = E\ = 
E\\ + E\ 2 , we can define the achievable rate-energy region as: 

C R _ E (P)±^(R,E) : R < logdet(I M + Hf 2 R^H 22 Q 2 ), 

E< I?* tr (Hy QjHg ) , tr (Qj ) < P, Q, y 0, j = lA (17) 

Here, because EH and ID operations in the IFC interact with each other, the boundary of the rate-energy 
region is not easily characterized and is so far unknown. The following lemma gives a useful insight into 
the derivation of the optimal boundary. 

Lemma 1: For Hu and H21, there always exists an invertible matrix T G C such that 

UgHnT = Eg 

VgH 21 T = I M , (18) 

where Ug* and Vg- are unitary and Sc is a diagonal matrix with ac,i > 0"g,2 o~g,m > 0. 

Proof: Because H 2 i has a full rank, by utilizing the generalized singular value decomposition [22], 
[25], we can obtain an invertible matrix T' such that 

UgHnT' = S A 

V£H 21 T' = S B , 
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where Ug and Vg are unitary and and are diagonal matrices with 1 > ua,i > ca,2 >,-••>> 
ca,m > and with < o"b,i < 05,2 <, < 0\b,m < 1, respectively. Here, i + i = 1. Therefore, 
by setting T = T'S^ 1 , we can obtain £[8]> with £ G = E^E^ 1 . ■ 
Without loss of generality, we set 

Qi = TXX H T H , (19) 

where X € C Mxm has the SVD as 

X = u^vf 

with T, x = diag{a Xtl , ...,a X!rn } and a X;1 >, > a x ^ m . Here, 

i=l 

where P' is a normalization constant such that ir(TXX^T^) < P is satisfied. We then have the 
following proposition. 

Proposition 2: In the high SNR regime, the optimal Qi at the boundary of the achievable rate-energy 
region has a rank one at most. That is, ranfe(Qi) < 1. 

Proof: First, let us consider the boundary point (R, E) of the achievable rate-energy, in which 
E < tr(Hi2Q2H^). Then, because the first transmitter do not need to transmit any signals causing the 
interference to the ID receiver (the second receiver), Qi = is optimal. That is, ranfc(Qi) = 0. 

For E > tr (H^C^Hfr,), let there be Qi with m = rank(Qi) > 1 which corresponds to the boundary 
point (R, E) of the achievable rate-energy. Then, given the harvested energy E (the boundary point) and 
Q2, the covariance matrix Qi exhibits 

£ = logdet(I M + Hf 2 R:iH 22 Q 2 ) (21) 

with 

ir(HnQiHg) = ^u, (22) 

where E\\ = E — tr '(H^G^H^). Because of Sylvester's determinant theorem [26] (det(I + AB) = 
det(I + BA) ), by substituting R_ 2 = 1m + H^iQiH^ into (|2~TT ). we can rewrite dill as 

R = logdet^ + HaaQaHg^M + HaiQxH^)" 1 ) 

= logdet((I M + H 2 iQiHg)(I M + HaxQiH^)- 1 + H 22 Q 2 Hf 2 (I M + HaiQiHg)" 1 ) 

= logdet((I M + HaiQiHg)" 1 + (H 2 iQiHg + H 22 Q 2 Hf 2 )(I M + H^QiH^)" 1 ) 

= logdet(I M + H 2 iQiHg + H 22 Q 2 Hf 2 ) - logdet(I M + H 2 iQiHg). (23) 
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Let us define 7n 2 = rank(Q2) and consider 771,2 > fn without loss of generality. From Lemma 1 and 
(fl9l ), d23l and (1221 can be respectively rewritten as 



R = logdet(I M + V G XX H Vg + H 22 Q 2 H|) - logdet(I M + V G XX fl Vf), 

= logdet(I M + XX" + Vf H 22 Q 2 Hf 2 V G ) - logdet(I M + XX H ), (24) 



and 



^(HuQiHfO = tr(HiiTXX H T H Hf x ) = tr(XJ G V G XX H i: H XJ%) 

m M 

= tr(S G XX^S G ) = ^a^.(^a Gi >F)| 2 ) = J E 11 , (25) 

j=l i=l 

where ui is the (i, j)th element of XJ X . From the interlacing theorem (Theorem 3.1 in [27]), (l24l can 
be further rewritten as 

(m m 2 \ m 

n(i+< 4 +N 2 ) n (1 + «j) 1 - io g + o^,*), 
i=l j=m+l I i=l 

(m m 2 \ m 

n«i+^ 2 ) n -iogn(i+4). ^ 
i=l j=m+l / i=l 

where re| is the interlaced value due to H22Q2H22. That is, <7y jma < ft| < cr^j, j = 1, m 2 , where <t 2 x 
and o" 2 m2 are the largest and the smallest eigenvalues of H22Q2H22. Note that the last approximation 
in (l26l ) is from the high SNR regime (i.e., large power P such that log(l + P) « log(P)), where cr^j 
for all i = 1, ...,7772 are linearly proportional to P resulting in 

Kj cc P for j = 1, .., 77i 2 . (27) 

Because <jg,i ctg,m > and 

M 
i=l 
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if there exists m > 1 such that (|25T ) is satisfied with d20b . we can find Q[ with rank(Q' 1 ) = 1 satisfying 
(l25l) . In addition, (l26l ) can be rewritten as: 

J II II ■ ill' ■ I i ! 

P 



(28) 




(29) 



FT"' /li ■ I ' (30) 

rii=i( i +<i) / 

The approximation in (|29l is from d27l l with a large P. That is because a^. i is negligible with respect 
to nf when E is finite. From (l20l . n^i(l + a \ i) i n tne denominator of (l30l has the minimum value 
when m = 1. In other words, if Qi with m > 1 exhibits (P, E), then we can find Q[ with m = 1 such 
that (P', E) with ^' > R in the high SNR regime, which contradicts that the point (R, E) is a boundary 
point. ■ 
Remark 1: Note that when the required harvested energy E (more precisely, En) is large, both a\ i 
and (Ty i are linearly proportional to P resulting in 



Then, (1281 ) becomes: 



oc P for i = 1, ...,m,j = 1, ..,m 2 . (31) 



R - log 1 + ^ n A • 02) 




Therefore, 



P oc log P m i- m = ( m2 - m ) bg P, (33) 



which implies that in the high SNR regime with large harvesting energy E, the achievable rate is linearly 
proportional to (rri2 — m). Then, we can easily find that it is maximized when m = 1. Note that it can 
be interpreted as the degree of freedom (DOF) in the IFC [28], in which by reducing the rank of the 
transmit signal at the first transmitter, the DOF at the second transceiver can be increased. 

Remark 2: Intuitively, from the power transfer point of view, Qi should be as close to the dominant 
eigenvector of HfjHn as possible, which implies that the rank one is optimal for power transfer. From 
the information transfer point of view, when SNR goes to infinity, the rate maximization is equivalent to 
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the DOF maximization. That is, a larger rank for Qi means that more dimensions at the second receiver 
will be interfered. Therefore, a rank one for Qi is optimal for both information and power transfer. 

When each node has a single antenna (M = 1), the scalar weight at the jth transmitter can be written 
as ■sJPjeP j or simply, Qj = Pj. The achievable rate-energy region can then be given as 

C R . E {P)^ UR,E) : R < log(l + xgig^), 

E<P x \h lx \ 2 + P 2 \h l2 \ 2 ,Pj<P,j = l,T>. (34) 

From (l34l . we can easily find that P2 = P at the boundary of the achievable rate-energy region. That is, 
the second transmitter always transmits its signal with full power P. Therefore, the optimal transmission 
strategy for M = 1 boils down to the power allocation problem of the first transmitter in the IFC. 

From Proposition [2] when transferring the energy in the IFC, the transmitter's optimal strategy is 
either a rank-one beamforming or no transmission according to the energy transferred from the other 
transmitter, which increases the harvested energy at the corresponding EH receiver and simultaneously 
reduce the interference at the other ID receiver. Even though the identification of the optimal achievable 
R-E boundary is an open problem, it can be found that the first transmitter will opt for a rank-one 
beamforming scheme. Therefore, in what follows, we first design two different rank-one beamforming 
schemes for the first transmitter and identify the achievable rate-energy trade-off curves for the two-user 
MIMO IFC where the rank-one beamforming schemes are exploited. 

B. Rank-one Beamforming Design 

1) Maximum-energy beamforming (MEB): Because the first receiver operates as an energy harvester, 
the first transmitter may steer its signal to maximize the energy transferred to the first receiver, resulting 
in a considerable interference to the second receiver operating as an information decoder. 

From Proposition |2j the corresponding transmit covariance matrix Qi is then given by 

Qi =Pi[V 11 ] 1 [V 11 ]f, (35) 

where Vji is a M x M unitary matrix obtained from the SVD of Hn and < P\ < P. That is, 
Hn = UiiEiiVf^, where En = diag{an : i, cr U)M } with 011,1 > ••• > 0"n,M- Here, the energy 
harvested from the first transmitter is given by Picr^ 1 . 
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2) Minimum-leakage beamforming (MLB): From an ID perspective at the second receiver, the first 
transmitter should steer its signal to minimize the interference power to the second receiver. That is, from 
Proposition [2j the corresponding transmit covariance matrix Qi is then given by 

Qi=Pi[V 21 ] M [V 2 i]^, (36) 

where V21 is a M x M unitary matrix obtained from the SVD of H21 and < P\ < P. That is, 
H 2 i = U2iS 2 iV^, where S 2 i = diag{a 2 i tl , cr 2 i,Af} with 0-21,1 > ... > o- 2 i,m- Then, the energy 
harvested from the first transmitter is given by P\ ||Hn [V 2 i]m|| 2 - 

C. Achievable R-E region 

Given Qi as in either (|33T l or (|36l , the achievable rate-energy region is then given as: 

Cr-e(P) = UR,E):R = R 2 ,E = E u + E 12 , 

R 2 < log det(I M + H^R:^H 22 Q2), E12 < tr (H 12 Q 2 Hf 2 ) , 

ir(Q 2 ) <P,Q 2 ^o|, (37) 

where 

{ P x o- 2 n 1 for MEB 
E n = { 11,1 , (38) 

{ Pi||Hii[V 2 i]m|| 2 for MLB 

and 

\ iM + AH^VnMVnlfHg for MEB 

k_ 2 = < rr (3y) 

[ Im + Pi^!i, m [ U 2i]m[U 2 i]^ for MLB 

Note that because a ul > ||Hn [V 2 i]m|| 2 , the energy harvested by the first receiver from the first 
transmitter with MEB is generally larger than that with MLB. 
Due to Sylvester's determinant theorem, R 2 can be derived as: 

R 2 = logdet(I M + Hf 2 Rl^H 22 Q 2 ) 

= logdet(lM + R: 2 /2 H 2 2Q2H|R:2 /2 ). (40) 

— 1/9 

Accordingly, by letting H 22 = R_ 2 H 22 , we have the following optimization problem for the rate-energy 
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region of (I37H 4 

(P3) maximize log det(I M + H 22 Q 2 Hg) (41) 
Q2 

subject to tr(Hi 2 Q 2 Hg) > max(P-Pn,0) (42) 

tr(Q 2 ) < P, Q 2 >= 0, (43) 

where P can take any value less than P max denoting the maximum energy transferred from both 
transmitters. Here, it can be easily derived that E mgx is given as 



E 

"max 



^11,1 + ^2,1) f° r MEB 
^(||Hn[V2i]jif|| 2 + a? 2il ) for MLB 



where o\ 2 x is the largest singular value of Hi 2 and it is achieved when the second transmitter also steers 
its signal such that its beamforming energy is maximized on the cross-link channel Hi 2 . That is, 

Q 2 = P[Vi 2 ]i[V 12 ]f. (45) 



Note that the corresponding transmit signal is given by x 2 (n) = vP[V 12 ]i,s 2 (n), where s 2 (n) is a random 
signal with zero mean and unit variance. Therefore, when s 2 (n) is Gaussian randomly distributed with 
zero mean and unit variance, which can be realized by using a Gaussian random code [30], the achievable 
rate is given by R 2 = logdet(I M + PH 22 [Vi 2 ]i [Vi 2 ]f Hg). 

Note that because En in (|42~1) and H 22 in (|4TT) are dependent on Pi (< P), we identify the achievable 
R-E region iteratively as: 

Algorithm 2. Identification of the achievable R-E region: 



1) Initialize n = 0, P 1 {0) = P, 



(0) p i °ni for MEB 

Pr\\H u {V 21 ] M \\ 2 for MLB 



and 



(0) / iM + PfWtVnHVulfH* for MEB 

lv, 2 = \ • (4 ' ) 

[ iM + PfVl^fU.^MtUs!]^ for MLB 

4 The dual problem of maximizing energy subject to rate constraint can be formulated, but the rate maximization problem is 
preferred because it can be solved using approaches similar to those in the rate maximization problems under various constraints 
[10], [19], [29]. 
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2) For n = : N, 
iptii 

If tr(H 12 Q^U^) + E[f >E 



Solve the optimization problem (P3) for Q 2 n ^ as a function of E^ and R^- 



E - ir(H 12 Q 2 n) Hg 



pf n+1 > = max I = v-^-v. izw ^Q j ( (48) 



where ft 



CTf a x for MEB 



llHntV^A/ll 2 for MLB 
Then, = min{P,p[ n+l) ) and update E$ +1) and R^ 2 +1) with P^ 4 "^ similarly to dHJ) and 

dUl). 

3) Finally, the boundary point of the achievable R-E region is given as (R,E) = (logdet(lM + 
H 22 Qf- +1 >H£), + tr(H 12 Qf- +1 )Hf 2 )). 

In (l48l ). if the total transferred energy (ir (Hj^Qj^Hf^) + is larger than the required harvested 

energy E, the first transmitter reduces the transmit power Pi to lower the interference to the ID receiver. 
In addition, if the energy harvested by the first receiver from the second transmitter (tr(Hi 2 Q 2 n ^H^)) 
is larger than E, the first transmitter does not transmit any signal. That is, rank(Q) = as claimed in 
the proof of Proposition [2] 

(n) 

To complete Algorithm 2, we now show how to solve the optimization problem (P3) for Q 2 in Step 

(n) (n) 

2 of Algorithm 2. The optimization problem (P3) with E X1 and R_ 2 can be tackled with two different 
approaches according to the value of E, i.e., < E < En and En < E < E max . Note that we have 
dropped the superscript of the iteration index (n) for notation simplicity. For < E < En, (P3) becomes 
the conventional rate maximization problem for single-user effective MIMO channel (i.e., H 22 ) whose 
solution is given as 

Q 2 = WF(H 22 ,I M ,P), (49) 

resulting in the maximum achievable rate for the given rank-one strategy Qi. Here, the operator WF() 
is defined in ( fTOl ). 

For .En < E < -E'max> the optimization problem (P3) can be solved by a "water-filling-like" approach 
similar to the one appeared in the joint wireless information and energy transmission optimization with 
a single transmitter [10]. That is, the Lagrangian function of (P3) can be written as 

L(Q 2 ,A,m) = logdet(I M + H 22 Q 2 Hf 2 ) 

+A(ir(H 12 Q 2 Hf 2 ) - {E-Ex)) - (i(tr{Q 2 ) - P), 
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and the corresponding dual function is then given by [10], [29] 

g(X, n) = maxL(Q 2 , A, (i). (50) 

Here the optimal solution //, A', and Q2 can be found through the iteration of the following steps [29] 

1) The maximization of L(Q 2 , X,pL) over Q2 for given A, //. 

2) The minimization of g(X, //) over A, /i for given Q2. 

Note that, for given A, [i, the maximization of L(Q 2 , A,y«) can be simplified as 

maxL(Q2,A,/u) = logdet(I M + H 22 Q 2 Hg) - ir(AQ 2 ), (51) 
Q2CZ0 

where A = /xIm — AH^H^. Note that ( f5TT > is the point-to-point MIMO capacity optimization with a 
single weighted power constraint and the solution is then given by [10], [29] 

Q 2 = A-^V> 22 A'V'£A-^, (52) 

where V 22 is obtained from the SVD of the matrix H22 A" 1 / 2 , i.e., H22A- 1 / 2 = U 22 £ 22 V 2 f . Here, 

£ 22 = diag{a' 2 2 1, ^'22 m} w i tn ^'221 — ■■■ — ^'22 M — ^ anc ^ ^ = diag{p\, ...,Pm} with = 
(1 — 1/(T 22 i = 1, M. The parameters /i and A minimizing gr(A, /x) in Step 2 can be solved by the 
subgradient-based method [10], [31], where the the subgradient of g(X, fj,) is given by (tr (H^C^H^) — 

(E-E 1 ),P-tr(Q 2 )). 

V. Energy-regularized SLER-maximizing beamforming 

In Section IIV-B[ two rank-one beamforming strategies are developed according to different aims - 
either maximizing transferred energy to EH or minimizing interference (or, leakage) to ID. Note that 
in [22], [23], the maximization of the ratio of the desired signal power to leakage of the desired signal 
on other users plus noise measured at the transmitter, i.e., SLNR maximization, has been utilized in the 
beamforming design in the multi-user MIMO system. Similarly, in this section, to maximize transferred 
energy to EH and simultaneously minimize the leakage to ID, we define a new performance metric, 
signal-to-leakage-and-harvested energy ratio (SLER) as 

I ITT v ||2 

SLER = " n - " — — . (53) 

||H 2 iv|| 2 + max(£-Pi||Hn|| 2 ,0) 

Note that the noise power contributes to the denominator of SLNR in the beamforming design [22], 

[23] because the noise at the receiver affects the detection performance degradation for information 

transfer. That is, the noise power should be considered in the computation of beamforming weights. In 

contrast, the contribution of the minimum required harvested energy is added in SLER of (l53l) . because 
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the required harvested energy minus the energy directly harvested from the first transmitter is a main 
performance barrier of the EH receiver. Therefore, in the energy beamforming, the required harvested 
energy is considered in the computation of the beamforming weights^ Then, the SLER of d53l ) can be 
rewritten as 

v^HgHuv 



" (Ug-li 21 +max(E/P 1 - ||H n || 2 , 0)I M ) v ' 
The beamforming vector v that maximizes SLER of (1541 ) is then given by 



(54) 



(55) 

|| V|| 

where v is the generalized eigenvector associated with the largest generalized eigenvalue of the matrix 
pair (Hf^Hn, H21H21 + max(E/P\ — ||Hn|| 2 , 0)Im). Here, v can be efficiently computed by using a 
generalized singular value decomposition (GSVD) algorithm [23], [32], which is briefly summarized in 
Algorithm 3. 

Algorithm 3. SLER maximizing GSVD-based beamforming: 
1) Set K = 



H n 
H21 



g £3MxM 

max {E/Pi - ||Hu|| 2 ,0)Im 

2) Compute QR decomposition (QRD) of K (= [P^P^R), where [P^P^] is unitary and R 6 
C MxM is upper Q-ianguiar Here Pq e c 2MxM . 

3) Compute V a from the SVD of P a , i.e., (P a ) 1;M V a = S a . 

4) v = R _1 [V a ]i and then, v 



w 



Here, because 



K: 



as in [32], for Pi||Hn|| 2 < E 



H n 

H21 



max (E/Pi - ||Hii|| 2 ,0)l 



M 



[P«;P/3]R 



R 



max 



(£/Pi-||Hn|| 2 ,0) 



(56) 



H^vp+la^fiinH^vip.o) - However, for computational simplicity, 

2 r> iixt l|2 



Strictly speaking, the SLER can be defined as SLER 
the lower bound on the required harvested energy is added in the denominator of SLER from the fact that ||Hn v||^ < Pi ||Hn 
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which avoids a matrix inversion in Step 4 of Algorithm 3. Because Algorithm 3 requires one QRD of an 
3M x M matrix (Step 2), one SVD of an M x M matrix (Step3), and one (M x M, M x 1) matrix-vector 
multiplication (Step 4 with d56ll), it has a slightly more computational complexity compared to the MEB 
and the MLB in Section IIV-BI that need one SVD of an M x M matrix. 

Once the beamforming vector is given as (l55t . we can obtain the R-E tradeoff curve for SLER 
maximization beamforming by taking the approach described in Section IIV-CI Interestingly, from (|54l , 
when the required harvested energy at the EH receiver is large, the matrix in the denominator of (l53l 
approaches an identity matrix multiplied by a scalar. Accordingly, the SLER maximizing beamforming is 
equivalent with the MEB in Section HV-B.li That is, v becomes V^PilVii]!. In contrast, as the required 
harvested energy becomes smaller, v is steered such that less interference is leaked into the ID receiver 
to reduce the denominator of (l53l . That is, v approaches the MLB weight vector in Section IIV-B.2I 
Therefore, the proposed SLER maximizing beamforming weighs up both metrics - energy maximization 
to EH and leakage minimization to ID. 

Note that the SLER value indicates how suitable a receiving mode, {EH\, ID2) or (ID\, EH2), is 
to the current channel. This motivates us to propose a mode scheduling between (EH\, ID2) and (ID±, 
EH2). That is, higher SLER implies that the transmitter can transfer more energy to its associated EH 
receiver incurring less interference to the ID receiver. Based on this observation, our scheduling process 
can start with evaluating for a given interference channel and P, 



2 



SLEBS^ = max- H Hli yH 

v ||H 2 iv|| 2 + maa;(£;-P||Hii|| 2 ,0) 



(57) 



and 



(58) 



SLERW = max- 

v Hi 2 vF + max{E — .P||H22 1| ,0) 

If SLERW > SLERP), (EHi, ID 2 ) is selected. Otherwise, (ID X , EH 2 ) is selected. 

VI. Discussion 

A. The rank-one optimality in the low SNR regime for one ID receiver and one EH receiver 

Even though we have assumed the high SNR regime throughout the paper, in some applications such as 
wireless ad-hoc sensor networks, low power transmissions are also considered. The following proposition 
establishes the rank-one optimality in the low SNR regime. 

Proposition 3: Considering (EH\, ID2) without loss of generality, in the low SNR regime, the optimal 
Qi at the boundary of the achievable rate-energy region has a rank one. 
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Proof: Similarly to (1231 ). the achievable rate at the ID 2 receiver is given by 

R = logdet(I M + HsxQxHfi + H 22 Q 2 Hg) - logdet(I M + H 2 iQiHfi). (59) 

For a Hermitian matrix A with eigenvalues in (—1, 1), logdet(I + A) can be extended as [33] 

log det(I + A) = tr(A) - -tr(A 2 ) + -tr(A 3 ) + .... (60) 

2 3 

Because HjjQjH^ is Hermitian and positive definite, and its maximum eigenvalue is upper-bounded 
as A max (HjjQjH^) < A max (Qj)A max (HjjH^) [33], for sufficiently low transmission power, their 
maximum eigenvalues lie in (—1, 1). Accordingly, in the low SNR regime, R can be approximated as 

R « 4r(H 2 iQiHg + H 22 Q 2 Hf 2 ) - ir(H 2 iQiHg) 

= tr(H 22 Q 2 H|). (61) 

That is, the achievable rate is independent of the interference from the first transmitter (noise-limited 
system). Then, Qi at the first transmitter can be designed to maximize the harvested energy. Therefore, 
the optimal Qi at the boundary of the achievable rate-energy region is given by 

Qi = argmax ir(HnQHn). (62) 
Q 

Note that fr(HnQHn) < Po\ x x , where the equality is satisfied when Q = P[Vn]i[Vn]f as in (|35T ). 
Therefore, the optimal Qi at the boundary of the achievable rate-energy region has a rank one. ■ 
Note that, from (loTl ). Q 2 maximizing R is designed as 

Q 2 = argmax ir(H 22 QHg) (63) 
Q 

and the corresponding Q 2 is given by Q 2 = P[V 22 ]i [V 22 ]^, where V 22 is the right singular matrix of 
H 22 . That is, at the low SNR, the optimal information transfer strategy in the joint information and energy 
transfer system is also a rank-one beamforming, which is consistent with the result in the information 
transfer system [34], where the region that the beamforming is optimal becomes broader as the SNR 
decreases. 

B. Asymptotic behavior for a large M 

Note that Proposition [2] gives us an insight on the joint information and energy transfer with a large 
number of antennas describing a promising future wireless communication structure such as a massive 
MIMO system [35]-[37]. 
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Given the normalized channel H = ^yrH, where the elements of H are i.i.d. zero-mean complex 

l H l 

Gaussian random variables (RVs) with a unit variance, and Q = Pw H with a finite P and ||v|| = 1, 
we define R = 1m + HQH^ which is analogous to R_2 in ((2Tb of the proof of Proposition |2] Then, 
because det(R) = 1 + Fv^H^Hv and 

H^H ps ^gH^H rj I m , (Central limit theorem in [38]) 

when M goes to infinity, det(R) « 1 + P and it is independent from the beamforming vector v. 
Analogously, when M increases, the design of vi in Qi = Piviv^ at the first transmitter is independent 
from det(R_2) (accordingly, independent from H21). Therefore, when nodes have a large number of 
antennas, the transmit signal for energy transfer can be designed by caring about its own link, not caring 
about the interference link to the ID receiver. That is, for a large M, MEB with a power control becomes 
optimal because it maximizes the energy transferred to its own link. 

Remark 3: Interestingly, from Section I VI- A I and IVI-BI we note that, when the SNR decreases or the 
number of antennas increases, the energy transfer strategy in the MIMO IFC would be designed by only 
caring about its own link to the EH receiver, not by considering the interference or leakage through 
the other link to the ID receiver. In addition, massive MIMO effect makes the joint information and 
energy transfer in the MIMO IFC naturally split into disjoint information and energy transfer in two 
non-interfering links. 

VII. Simulation Results 

Computer simulations have been performed to evaluate the R-E tradeoff of various transmission 
strategies in the two-user MIMO IFC. In the simulations, the normalized channel Hy is generated 
such as Hjj = p^^ Hjj, where the elements of are independent and identically distributed (i.i.d.) 
zero-mean complex Gaussian random variables (RVs) with a unit variance. The maximum transmit power 
is set as P = 50W, unless otherwise stated. 

Table U lists the achievable rate and energy, (R, E), for single modes - (EH\, EH2) and (ID\, ID2). 
The harvested energy of {EH\, EH2) for M = 4 is larger than that for M = 2 and furthermore, the 
achievable rate of (ID\, ID2) for M = 4 is higher than that for M = 2. Note that the achievable rate 
of (EHi, EH2) and the harvested energy of (ID\, ID2) are zero. 

Fig. |2] shows R-E tradeoff curves for the MEB and the MLB described in Section ITV-BI when M = 4, 
an = 1 for i = 1, 2, and ay = 0.8 for i ^ j. The first transmitter takes a rank-one beamforming, either 
MEB or MLB, and the second transmitter designs its transmit signal as d45l ), d49l ), and ( [521 , described in 
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TABLE I 

The achievable rate and energy, (R, E), for single modes when M g {2, 4} 



Mode 


M = 2 


M = 4 


(EH U EH 2 ) 


(0, 262.98) 


(0, 359.57) 


{ID X , ID 2 ) 


(9.67, 0) 


(16.08, 0) 



Section IIV-CI As expected, the MEB strategy raises the harvested energy at the EH receiver, while the 
MLB increases the achievable rate at the ID receiver. Interestingly, in the regions where the energy is less 
than a certain threshold around 45 Joule/sec, the first transmitter does not transmit any signals to reduce 
the interference to the second ID receiver. That is, the energy transferred from the second transmitter is 
sufficient to satisfy the energy constraint at the EH receiver. 

The dashed lines indicate the R-E curves of the time-sharing of the full-power rank-one beamforming 
(either MEB or MLB) and the no transmission at the first transmitter. Here, the second transmitter 
switches between the beamforming on H12 as d45l ) and the water-filling as (|49l i in the corresponding 
time slots. For MLB, "water-filling-like" approach (l52l exhibits higher R-E performance than the time- 
sharing scheme. However, for MEB, when the energy is less than 120 Joule/sec, the time-sharing exhibits 
better performance than the approach (|52l . That is, because the MEB causes large interference to the ID 
receiver, it is desirable that, for the low required harvested energy, the first transmitter turns off its power 
in the time slots where the second transmitter is assigned to exploit the water-filling method as (|49l i. 
Instead, in the remaining time slots, the first transmitter opts for a MEB with full power and the second 
transmitter transfers its information to the ID receiver by steering its beam on EH receiver's channel 
H12 as (|45T ) to help the EH operation. In Fig. [3l we have additionally included the R-E tradeoff curves 
for MEB with rank(Qi) = 2 when the simulation parameters are the same as those of Fig. |2] Here, 
we can find that the MEB with rank(Qi) = 1 has superior R-E boundary points compared to that with 
rank(Qi) = 2. That is, even though we have not identified the exact optimal R-E boundary, for a given 
beamforming (MEB in Fig. 3), we can find that the beamforming with rank(Qi) = 1 has superior R-E 
boundary points compared to that with rank (Qi) = 2. 

In Fig. @1 we plot R-E tradeoff curves for M = 4 and P = 0.1. As observed in Section IVI-A1 at 
the low SNR, the MEB exhibits higher harvested energy than the MLB without any degradations in the 
achievable rate. Fig. [5] shows R-E tradeoff curves for M = 15. Compared to M € 4 (Fig. [2), the gap 
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Fig. 2. R-E tradeoff curves for MEB and MLB when M = 4, ce H = 1 for i = 1, 2, and ay = 0.8 for i j. 



between the achievable rates of MEB and MLB is relatively less apparent. As pointed out in Remark [3] 
of Section EH for low SNRs or large numbers of antennas in the MIMO IFC, the energy transfer strategy 
of maximizing the transferred energy on its own link exhibits wider R-E region than that of minimizing 
the interference to the other ID receiver. 

Fig. [6] shows R-E tradeoff curves for MEB, MLB, SLNR maximizing beamforming, and SLER 
maximizing beamforming when M = 4, am = 1 for i = 1,2, and aij = 0.8 for i ^ j. The R-E 
region of the proposed SLER maximizing beamforming covers most of those of both MEB and MLB, 
while the SLNR beamforming does not cover the region for MEB. Fig. [7] shows R-E tradeoff curves for 
an asymmetric case M t = 3 and M r = 4. We can find a similar trend with M t = 4 = M r = 4, but the 
overall harvested energy with M t = 3 and M r = 4 is slightly less than that with M t = 4 = M r = 4. 

Fig. [8] shows the R-E tradeoff curves for SLER maximizing beamforming with/without SLER-based 
scheduling described in Section [V] when (a) = 0.7 and (b) = 1 for i ^ j. Here, we set an = 1 
for i = 1, 2 and M = 2. Note that the case with a>ij = 0.7 has weaker cross-link channel (inducing less 
interference) than that with = 1. The SLER-based scheduling extends the achievable R-E region for 



May 1, 2013 



DRAFT 



24 



R-E tradeoff curves when M =M =4 
r t 

250 1 1 1 1 1 ■ — 




8 9 10 11 12 13 

Rate (bits/s/Hz) 

Fig. 3. R-E tradeoff curves for MEBs (ranfc(Qi) = {1, 2}) and MLB (ranfc(Qi) = 1) when M = 4, an = 1 for i = 1, 2, 
and djj = 0.8 for j 7^ j. 

both ajj 6 {0.7, 1}, but the improvement for atij = 1 is slightly more apparent. That is, the SLER-based 
scheduling becomes more effective when strong interference exists in the system. Note that the case with 
ctij = 1 exhibits slightly lower achievable rate than that with ctij = 0.7, while achieving larger harvested 
energy. That is, the strong interference degrades the information decoding performance but it can be 
effectively utilized in the energy-harvesting. 

VIII. Conclusion 

In this paper, we have investigated the joint wireless information and energy transfer in two-user 
MIMO IFC. Based on Rx mode, we have different transmission strategies. For single-operation modes 
- (ID\, ID2) and {EH\, EH2), the iterative water-filling and the energy-maximizing beamforming on 
both receivers can be adopted to maximize the information bit rate and the harvested energy, respectively. 
For (EHi, ID2), and (ID\, EH2), we have found a necessary condition of the optimal transmission 
strategy that one of transmitters should take a rank-one beamforming with a power control. Accordingly, 
for two transmission strategies that satisfy the necessary condition - MEB and MLB, we have identified 
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R-E tradeoff curves when M_=M ( =4 and P=0.1 
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Fig. 4. R-E tradeoff curves for MEB and MLB when M = 4 and P = 0.1. 

their achievable R-E tradeoff regions, where the MEB (MLB) exhibits larger harvested energy (achievable 
rate). We have also found that when the SNR decreases or the number of antennas increases, the joint 
information and energy transfer in the MIMO IFC can be naturally split into disjoint information and 
energy transfer in two non-interfering links. Finally, we have proposed a new transmission strategy sat- 
isfying the necessary condition - signal-to-leakage-and-energy ratio (SLER) maximization beamforming 
which shows wider R-E region than the conventional transmission methods. That is, we have found that 
even though the interference degrades the ID performance in the two-user MIMO IFC, the proposed 
SLER maximization beamforming scheme effectively utilizes it in the EH without compromising ID 
performance. 

Note that, motivated from the rank-one beamforming optimality, the identification of the optimal R-E 
boundary will be a challenging future work. Furthermore, the partial CSI or erroneous channel information 
degrades the achievable rate and the harvested energy at the receivers, which drives us to develop a robust 
rank-one beamforming. 



■ Max. energy beamforming 

■ Min. leakage beamforming 
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Max. energy beamforming 
Min. leakage beamforming 
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Fig. 5. R-E tradeoff curves for MEB and MLB when M = 15. 
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